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Abstract 

The process e"’'e“ —>■ "fXcj {J=l, 2) is studied via initial state radiation using 980 fb“^ of 
data at and around the T{nS) (n=l, 2, 3, 4, 5) resonances collected with the Belle detector at 
the KEKB asymmetric-energy e'*"e“ collider. No significant signal is observed except from b(2«S') 
decays. Upper limits on the cross sections between ^/s = 3.80 and 5.56 GeV are determined at 
the 90% credibility level, which range from few pb to a few tens of pb. We also set upper limits 
on the decay rate of the vector charmonium [■(/'(dOdO), "0(4160), and 0(4415)] and charmoniumlike 
[y(4260), y(4360), and y(4660)] states to xXcj- 

BAGS numbers: 14.40.Pq, 13.25.Gv, 13.66.Bc 
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In e^e~ annihilation, the energy region above the DD threshold is rich with vector 
charmonium and charmoniumlike states. Three charmoniumlike states with = 1 were 
discovered at B factories via initial state radiation (ISR) in the last decade: the R(4260) 
in e+e" —)■ 7r+7r“J/'0 [H, and the y(4360) and y(4660) in —)■ 7 r''' 7 r“' 0 ( 2 S') [ 1 , |^. 

Together with the conventional charmonium states '0(4040), '0(4160), and 0(4415), there 
are six vector states; the potential models predict only hve in this mass region [^. Some of 
these states show unusual properties that are inconsistent with charmonium . It is unlikely 
that all of these states are charmonia; some, perhaps, have exotic nature: a multiquark state, 
molecule, hybrid, or some other configuration. To improve our understanding of these states 
and the underlying QCD, it is important to investigate them using much larger data samples 
and new decay channels. 


For example, one can study radiative transitions between these states and lower char¬ 
monium states like the Xcj- The CLEO Collaboration used data taken during a scan of 
center-of-mass (CM) energies ^/s = 3.97 — 4.26 GeV to report upper limits on the cross 
sections of —)■ xXci and e+e“ —)■ xXc2 in three energy regions: the 0(4040) (-0s = 3.97- 
4.06 GeV), the 0(4160) (4.12-4.20 GeV), and y/s = 4.26 GeV . The limited statistics 
prevented them from measuring the line shape of e’''e“ —)■ xXcj- The BESIII experiment 
reports the upper limits on the cross sections of the reactions e’''e“ —)■ xXa and e+e“ —)■ 7 Xc 2 
at four energy points: -0s = 4.009, 4.230, 4.260, and 4.360 GeV With the full Belle data 
sample, we are able to study this process via ISR. 


In this paper, we report a study of the —)■ xXcj process using ISR events detected 

with the Belle detector at the KERB asymmetric-energy collider 10|. Here, XcJ 

is reconstructed in the 7<//0 hnal state and J /0 is reconstructed in the fjT final state 
alone (The background level is very high in the hnal state due to Bhabha events). The 
same hnal state 77 J/ 0 , has been previously analyzed at Belle and 0(4040) and 0(4160) 
were observed as 7 J/0 resonances [ll[. We study the full Belle dataset corresponding to 
an integrated luminosity of 980 fb“^. About 70% of the data were collected at the T(4S') 
resonance, and the remainder were taken at the other T(nS') (n=l, 2 , 3, or 5) states or at 
CM energies a few tens of MeV lower than the T (45*) or the T {nS) peaks. 

The event generator evtgen 1^ with the VECTORISR model is used to simulate the signal 
process e+e“ —>■ 7isrR 7isr7XcJ 7isr77T/0. The mass and width of V can be varied so 
that we can obtain the signal efficiency as a function of the vector meson mass. This model 
considers the leading-order (LO) quantum electrodynamics (QED) correction only and thus 
higher-order corrections should be estimated and properly taken into account. The dedicated 
ISR generator phokhara [l^ has the next-to-Ieading-order (NLO) QED correction but does 
not contain the mode of interest. However, the process —)■ 7 isrR —^ 7isr 7J/0 can be 

generated with phokhara and this allows us to estimate the NLO correction effect in the 
mode under study by comparing the results from the two generators in the analysis of the 
7 J/0 mode. All generated events are passed through the GEANT3 
simulation and then the standard reconstruction. 


3 based detector 


For a candidate event, we require two good charged tracks with zero net charge. The 
impact parameters of these tracks perpendicular to and along the beam direction with 
respect to the interaction point are required to be less than 0.5 cm and 5.0 cm, respectively. 
The transverse momentum of the leptons is required to be greater than 0.1 GeV/c. For 
each charged track, information from different detector subsystems is combined to form a 
likelihood Ci for each particle species (i) [3]- For muons from J/ip ^ one of the tracks 
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is required to have the muon identihcation likelihood ratio TZu = r Lr > 0.95; in addition, 
if one of the muon candidates has no muon identihcation (ID) information [ig, the polar 
angle of each muon candidate in the 7XcJ CM system is required to satisfy | cos 0^1 < 0.75. 
The lepton ID efficiency is about 87% for J/ijj —)■ 

A photon candidate is an electromagnetic calorimeter cluster with energy 77(7) > 50 MeV 
that does not match any charged tracks. The photon is labeled as the ISR photon when 

its energy in the e+e“ CM frame exceeds 3 GeV (corresponding to MIjXcj] < 7 GeV/c^, 

the maximum non-ISR photon energy being about 3 GeV) and this photon is excluded 
when reconstructing xXcj candidates. We also require at least two additional photons, 
each with energy in the laboratory frame greater than 0.25 GeV. Among these, we select 
the two with the highest energy in the laboratory system and denote these as jh and 7/ 
(with > E^^). The detection of the ISR photon is not required; instead, we require 
— 1 (GeV/c^)^ < < 2 (GeV/c^)^, where is the square of the mass recoiling against 

the 7Xcj system. The distribution of is shown in Fig. [T] 



FIG. 1. Missing mass squared distribution with M{'yi'yhJ/'4’) < 5.56 GeV/c^. 

Fig. [2] shows the invariant mass [Mfi~)] distribution for events that survive 

the selection criteria and with the XilhJ/iJ invariant mass u~) — 

^~) + mj/^] less than 5.56 GeV/c^, where mj/^ is the nominal mass of the J/'ip |l8j |. 
A pair is considered as a. J/pj candidate if M(u'^u~) is within ±45 MeV/c^ (the mass 
resolution being 15 MeV/c^) of the J/ij) nominal mass [181]. The J/ip mass sidebands are 
dehned as M(/i+/i“) G [3.172, 3.262] GeV/c^ or [2.932, 3.022] GeV/c^, which are twice as 
wide as the signal region. 

To reject the background from e'''e“ 7isR7(7r°)T/'0 events with 7 or decaying into 

two photons, we require that the invariant mass of the two photons, M{xi)i be outside the 
7] mass region of [0.50, 0.58] GeV/c^, the 7r° mass region and the low-invariant-mass region 
< 0.20 GeV/c^. Figure [3] shows the invariant mass distribution of M{''yJ/'ip) (with two 
entries per event for M{xhJ/pj) and for events with M{■jilhJ/i’) < 5.56 GeV/c^. 

Here, IJ'~) — fi~) + mj/^. We observe xa and Xc 2 signals but 

no evidence of Xco- We divide the XcJ niass region into [3.48, 3.535] GeV/c^ for Xci and 
[3.535, 3.58] GeV/c^ for Xci- 

Figured] shows the M{'ynhJ distribution after applying all the selection criteria above. 
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FIG. 2. Invariant mass distribution of . The shaded area in the middle is the J/V’ signal 
region, and the shaded regions on either side are the J/ij) mass sidebands. 



M(yJ/\(/) (Gev/c^) 

FIG. 3. Invariant mass distribution of l'i> for candidate events with M/if:) < 5.56 GeV/c^. 
The shaded histograms show the Xci ([3.48, 3.535] GeV/c^) and Xc 2 ([3.535, 3.58] GeV/c^) regions. 


We see a clear 'ip{2S) signal but no significant signal in the higher mass region. The clear Xcj 
and 'ip{2S) signals allow us to measure the product branching fractions B['ip{2S) —?■ xXcj] x 
^[Xcj —t tJ/V'] {J = 1, 2). By contrast, in the region G [3.80, 5.56] GeV/c^, we 

set an upper limit on the production cross section of XXcj- 


The potential backgrounds are also shown in Fig. IH Besides the non-J/ip background, 
which also appear in the J/ip mass sidebands, there are three additional backgrounds: 
7 isr-^/^, TiSKTr^Tr^J/V^, and XiSKVJ/i>- Of course, e+e“ -)■ 7 isrV'( 2S') with ip{2S) -)■ 


e' e 


IXcj will be a background in the analysis of the XilhJ/'ip high-mass region. The ISR J/ip and 
ip{2S) samples are generated according to the theoretical calculation of the production cross 
sections 17] with the world-average resonant parameters as input [l^. For the other modes, 
we use the cross sections of e+e“ rjj/ip [Tl| and —)■ J/ip [HI and assume that 


a 


{e^e — TT^TT^J/ip) = ^a(e~^e —>■ 7r’''7r J/ip). All these samples are generated using the 
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FIG. 4. Invariant mass distribution of "fnhJ/'4’- The background from the tail of the V'(2S') is 
plotted only for M/ijj) > 3.75 GeV/c^ and M{'yi^hJ/'ip) < 3.65 GeV/c^. The dots with error 
bars are data while the shaded histograms represent different sources of background modes. 


PHOKHARA generator [Isj and are normalized to the integrated luminosity of the full data 
sample. The background contribution practically saturates the mass spectrum above the 
'ijj{2S) peak. 

To measure the '4>{2S) —)■ yXcj branching fractions, we dehne the 'ip{2S) signal region as 
3.65 GeV/c^ < /iIj) < 3.72 GeV/c^. The distribution of the energy of the less ener¬ 
getic photon in the GM system is shown in Fig. [5l Glear signals due to xa and Xc 2 

are observed with very low background and we ht this photon energy distribution to extract 
the corresponding yields. The Xcj signal shapes are obtained from Monte Garlo simulated 
signal samples convolved with a corresponding smearing Gaussian function to compensate 
for the resolution difference between data and Monte Garlo simulation; the background is 
parameterized as a hrst-order Ghebyshev polynomial. The resulting ht function is shown in 
Fig. [5] and the ht yields 340 ± 20 Xci and 97 ± 12 Xc 2 signal events. 

From the world-average 4>{2S) resonant parameters [l8|, we calculate cr[e’''e“ —)■ 7 isr'^( 2S')] 
= (14.25 ± 0.26) pb [T^ and thus expect 13.9 x 10® ISR produced '^(25') events in the full 
Belle data sample of 980 fb“h With the efhciencies of 1.4% and 0.7% for the Xci and 
Xc 2 modes, respectively, from the MG simulation, we obtain B[il>{2S) —?■ xXci] x ^{Xci 
xJ/ip) = (2.92 ± 0.19)% and B['ijj{2S) —)■ xXc 2 \ x B{xc 2 = (1-65 ± 0.21)%. Here, 

the errors are statistical only. These results are consistent with the PDG values [l8|. 

The M{xi'jhJ /distributions above the '4>{2S) signal region for xXci and xXc 2 candidate 
events as well as their sum are shown in Fig. [61 together with the background estimation 
from the J/if) mass sidebands and the MG simulated background modes with a genuine 
J/-^. No signihcant signal is observed in either the xXa or xXc 2 mode. As the background 
estimation is limited to the known channels, it only serves as a lower limit of the true 
background. In calculating the upper limits of the xXcJ production cross section, we consider 
the estimated-background events from the observed signal candidates. This results in a 
conservative estimate of the upper limit of the signal and hence a conservative estimate for 
the cross section. 

There is cross contamination between the Xci and Xc 2 signals due to the mass resolution. 
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(GeV) 


FIG. 5. Energy distributions of the low energy photon in the 'yi'yhJ/4’ CM system for events 
in the 'ip{2S) mass region. Dots with error bars are data and histograms are MC samples. The 
blue solid line is the best fit, the red dashed line is the shape of the total background determined 
from the fit, and the purple dot-dashed line is the MC signal shape convolved with a Gaussian 
function. The shaded histogram shows the total background as determined from J/'0 sidebands 
and simulations. 

as can be seen from Fig. [3l and this is taken into account as follows. The yields of observed 
Xci and Xc 2 events (denoted as and respectively) are expressed as 


( ^ohs ] = f £21 ] ( X B{xci iJ/'tp) X p+p )\ , ( ^bkg \ ^ 

V / V ^12 £22 ) V X ^iXc2 iJ/'ip) X B{J/iIj /i+p-) ) \ )' 

In these equations, eij {i, j = I, 2) is the efficiency of produced xa to be reconstructed 
in the Xcj signal region; and represent the total numbers of Xci and Xc2 events 
produced in data, respectively; B is the world-average branching fraction for the given 
process [l^; and n^kg and n^kg represent the numbers of non-Xcj background events for Xci 
and Xc 2 , respectively, which are the sum of the normalized J/'0 mass sideband background 
and the MC simulated J/-^, and 7 isr' 0(2S') background, as 

shown in Fig. | 6 l The efficiency curves cn and 622 , also shown in Fig. [ 6 l are not monotonic 
between 3.9 GeV/c^ < m{xXcj) < 4.2 GeV/c^. This is due to the fact that the energies of 
the two photons are almost the same in this mass region. 

We use the maximum likelihood method to determine upper limits on the numbers of 
produced xXcj events, and and thus on the upper limits of the production cross 
sections of e’''e“ —)■ xXcJ- The likelihood is constructed as follows. For each possible pair of 
the A^Xci ]\fXc 2 values, the numbers of the expected signal events, and are 

/ jyxci \ _ / eii €21 \ X B{xci iJ/i}) X B^J/iIj pV) \ 

Z/Xc2 j €22 ) \ X B{Xc2 ^ iJ/^) X B{J/Ip ^ ^ 

Taking into account the background contribution, the numbers of expected events in the 
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M(yj^J/(|/) (GeV/c^) 


FIG. 6. Invariant mass distributions of JXcJ candidates. Shown from top to bottom are xXcij 
7Xc2; and their sum. Dots with error bars are data, the shaded histograms are the simulated 
backgrounds and J/V’ sidebands, and the solid lines are the efficiency curves. 


signal regions, denoted as and for Xci and Xc 2 , respectively, are 





and the probability of observing 


Xcl 

obs 

Xc2 

obs 


events in data is 


( 3 ) 






n 


Xcl I 
obs ■ 


n 


Xc2 I 
obs ■ 


( 4 ) 


The nncertainty in the background estimation is considered by sampling in Eq. (E]). 
By fitting the normalized background distribution, the mean value and the uncertainty 
of the background level are obtained. The background yield n^^g is varied assuming it 
follows a Gaussian distribution with this mean value and the uncertainty as the standard 
deviation. The systematic error of the measurement, which corresponds to an uncertainty in 
the expected number of events, follows a Gaussian distribution with a mean value and 
a standard deviation x (Tsys, where Usys is the total relative systematic error (13.4%), 
described below. This is also considered by varying in Eq. (jl]). 

The summation of random-sampled , considering the uncertainty in back¬ 

ground estimation and the systematic errors, forms the final likelihood function 
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( 5 ) 


Here, N is the number of samplings, +^^bkg,/ "'^^ere 

^bkg,i) n are the numbers of events obtained from the corresponding Gaussian 

distributions. The subscript k represents the fc-th sampling for the expected number of Xci 
signal events z/^=b The other subscripts /, m and n have parallel meanings. By letting 
and run over all the possible values from 0 to inhnity independently, we obtain 
the likelihood in the plane. The likelihood L{N^‘^^) can be obtained from this 

two-dimensional likelihood function by integrating over the variable From this, we 

obtain the upper limit on at the 90% credibility level (C.L.) 0 and convert this into 
the upper limit on a{e'^e~ —)■ 7Xci)- The upper limit on a{e'^e~ 7Xc2) is determined in 

a similar manner. The hnal upper limits are shown in Fig. [7] and are around a few pb to a 
few tens of pb. We also show the CLEO and BESIII results in Fig. [7] for comparison. The 
measured upper limits are more stringent than the CLEO results at ^/s = 3.97 — 4.06 GeV 
and y/s = 4.26 GeV. The large data samples collected by BESIII at y/s = 4.009, 4.230, 
4.260, and 4.360 GeV provide stronger upper limits at these energy points. The values of 
the upper limits measured here are listed in Table [T1 

We extract the transition rate of the vector charmonium and charmoniumlike states to 
7XcJ by htting the distributions in Fig. [6l We use a Breit-Wigner function for the signal 
and a hrst- or second-order polynomial function for the background. While doing the ht, 
the mass and total width are hxed to the world average-values [isj and Fee x H(i? —)■ xXcj) 
is scanned from zero to a large number at which the probability is less than 1.0% of the 
largest value. Normalized probability density functions are derived from such a scan. These 
probability density functions then give the upper limits at 90% C.L. as listed in Table UTl 
Taking Fee['0(4O4O)] and Fee['^(4415)] from the world average-values [l^ and Fee['0(416O)] 
from the BES II measurement [1^, we set the upper limits on the branching fractions for 
these three conventional charmonium states as listed in Table HTTl Taking F^p \Y ('4260)] x 
B[Y (4260) J/tp] = (6.4 ± 0.8 ± 0.6) eV or (20.5 ± 1.4 ± 2.0) eV [l9[ (there are two 

solutions for the best ht in this mode, and there are also two solutions in the Y (4360) and 
Y (4660) cases below) ), Fee[V(4360)] x B[Y (4360) —?■ 7r+7r“'0(2S')] = (10.4± 1.7± 1.4) eV or 
(11.8T1.8T1.4) eV fl, and Fee[V(4660)] x H[V(4660) ^ 7r+7r-^/;(2^)] = (3.0T0.9T0.3) eV 
or (7.6 ± 1.8 ± 0.8) eV j^, we set the upper limits on the ratios of the branching fractions 
as shown in Table IIVI The mass and width of the vector charmonium and charmoniumlike 
states, the background shape, and the ht range are varied in the ht to estimate the systematic 
uncertainties. The largest upper limit from these tests is taken as the hnal result. The 
total uncertainties from the reference processes and the systematic errors are considered by 
assuming they are Gaussian errors. 

The following sources of systematic uncertainties are considered in the a{e'^e~ XXcj) 
upper-limit determination. The uncertainty in the tracking efficiency for tracks with angles 
and momenta characteristic of signal events is about 0.35% per track Ea and is additive. 
The uncertainty due to particle identihcation efficiency is 1.9%. The uncertainty oi J/ip mass 


^ In common high energy physics usage, this Bayesian interval has been reported as “confidence interval” 
which is a frequentist-statistics term. 
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TABLE I. Upper limits on the e'’'e 7Xcj cross sections. 


^/s (GeV) 

Xci (pb) 

Xc2 (pb) 

(GeV) 

Xci (pb) 

Xc2 (pb) 

3.80-3.84 

80 

134 

4.68-4.72 

8 

14 

3.84-3.88 

37 

90 

4.72-4.76 

8 

18 

3.88-3.92 

35 

no 

4.76-4.80 

11 

15 

3.92-3.96 

27 

40 

4.80-4.84 

9 

18 

3.96-4.00 

12 

21 

4.84-4.88 

15 

11 

4.00-4.04 

34 

53 

4.88-4.92 

11 

14 

4.04-4.08 

29 

45 

4.92-4.96 

10 

10 

4.08-4.12 

46 

54 

4.96-5.00 

4 

21 

4.12-4.16 

27 

53 

5.00-5.04 

8 

13 

4.16-4.20 

10 

63 

5.04-5.08 

13 

13 

4.20-4.24 

36 

35 

5.08-5.12 

11 

7 

4.24-4.28 

14 

17 

5.12-5.16 

9 

7 

4.28-4.32 

19 

38 

5.16-5.20 

5 

17 

4.32-4.36 

16 

20 

5.20-5.24 

14 

9 

4.36-4.40 

8 

22 

5.24-5.28 

7 

6 

4.40-4.44 

14 

34 

5.28-5.32 

6 

8 

4.44-4.48 

11 

22 

5.32-5.36 

4 

16 

4.48-4.52 

11 

21 

5.36-5.40 

6 

14 

4.52-4.56 

7 

12 

5.40-5.44 

4 

10 

4.56-4.60 

16 

13 

5.44-5.48 

8 

8 

4.60-4.64 

6 

26 

5.48-5.52 

8 

8 

4.64-4.68 

12 

20 

5.52-5.56 

4 
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TABLE 11. Upper limits on Tee x B at the 90% C.L. 



Xci (eV) 

Xc2 (eV) 

ree[V'(4040)] X ^[-0(4040) 7Xcj] 

2.9 

4.6 

ree[V'(4160)] X B[V’(4160) ^ 7Xcj] 

2.2 

6.1 

ree[V'(4415)] xB[V’(4415) ^7Xcj] 

0.47 

2.3 

ree[V(4260)] X B[V(4260) ^ 7Xcj] 

1.4 

4.0 

ree[V(4360)] X B[V(4360) ^ -fXcj] 

0.57 

1.9 

ree[V(4660)] X B[V(4660) ^ xXcj] 

0.45 

2.1 


TABLE III. Upper limits on branching fractions B{R —)• 'yXcj) at the 90% C.L. 


Resonance 

7Xci (10 ^) 

7Xc2 (10 ^) 

^/;(4040) 

3.4 

5.5 


6.1 

16.2 

V'(4415) 

0.83 

3.9 


11 
























































{s (GeV) 


FIG. 7. Measured upper limits on the e"’'e“ —>• 7Xcj cross sections at the 90% C.L. for Xci (top) 
and Xc 2 (bottom). The solid dots show the Belle measurements, the solid triangles are the resnlts 
from CLEO and the blue squares are from BESIII. 


TABLE IV. Upper limits on branching fraction ratios at the 90% C.L. The two upper limits 
correspond to the two solutions in the reference processes. 


Resonance 

IXcl 

lXc2 

Biy(4260)^7XcjJ 

0.3 or 0.07 

0.7 or 0.2 

B[y(4260)-i>7r+7r-J/i/)l 

BIK(4360)^7XcjJ 

0.06 or 0.05 

0.2 or 0.2 

B\Y (4360)— >-7r+7r“i/>(25)l 

Biy(4660)^7XciJ 

0.2 or 0.07 

0.9 or 0.3 

B[y(4660)->7r+7r-i/)(2S)l 
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and XcJ mass requirements are estimated using the 'ip{2S) sample in the same analysis and 
they are found to be 1% and 1.3%, respectively. The generator evtgen is used in generating 
signal MC events. In this generator, however, only one ISR photon is allowed and the higher- 
order ISR effect should be estimated and corrected. This effect is studied by using a control 
sample 7isr'0(2*S') with 'ip{2S) decaying into rjj/'il). This process can be generated 

with both EVTGEN and phokhara, a generator with higher-order ISR corrections. We 
assume that the correction factor obtained in this mode is the same as in the mode under 
study, and 9.0% is taken as the systematic error, corresponding to the uncertainty in the 
difference between the measured B{'ijj{2S) —)■ JXcj —t and the world average jisj. 

Taking the statistical error of the MC samples and the possible uncertainty in simulating the 
angular distributions of the full decay chain xXcj mto account, we quote a total 

uncertainty due to the generator as 12%. Belle measures luminosity with 1.4% precision 
and the trigger efficiency is about 91% with an uncertainty of 2%. Errors on the branching 
fractions of the intermediate states are taken from Ref. [l8| with a systematic error of 4.5%. 
Assuming that these systematic error sources are independent, the total systematic error is 
13.4%. The systematic uncertainty is considered in the upper limits shown in Tables HI— IIVI 

In summary, using the full Belle data sample, we measure the —)■ xXcJ process 

via initial state radiation. For the CM energy between 3.80 and 5.56 GeV, there are no 
signihcant —)■ xXa and xXc 2 signals. The upper limits on the XXcj production 

cross sections, which range from a few pb to a few tens of pb, are set for the first time and 
are listed in Table HI We also set upper limits on the decay rate of the vector charmonium 
and charmoniumlike states to xXcJ- This information may help in understanding the nature 
of these vector states. 
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